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1.  INTRODUCTION  • 

The  quadrupole  mass  filter  (QMF)  is  ubiquitous  in  military  and 
civilian  mass  spectrometry  laboratories.  This  instrument  is  used  by 
chemists  to  determine  the  mass  distribution  associated  with  a  sample.  As 
illustrated  in  figure  1,  a  QMF  has  four  pole  pieces  and  each  of  these  is 
carefully  shaped  to  follow  the  equation  of  a  hyperbola  i.e.  with  the  co¬ 
ordinate  system  shown  the  equations  of  the  pole  pieces  are  given  by 
(x2  -  y2)  =  ±r o.  In  a  typical  QMF,  rQ  and  l  are  .28  cm  and  15  cm 
respectively  (1).  Detailed  analysis  shows  (2)  that  mass  filtering  action 
is  achieved  by  imposing  a  voltage  difference  <j>  =  -  <f>  cos  tot 

between  adjacent  electrodes.  The  driving  frequency  V  =  (<jj/2tt)  is 

typically  2.5  MHz.  The  voltages  tj>  and  are  refered  to  as  the  DC 
and  AC  components  of  the  applied  voltage  respectively.  The  result  of 
detailed  analyses  (1,2)  is  that  the  mass  resolution  of  a  QMF  is  determined 
by  the  ratio  *  while  the  absolute  magnitude  of  the  voltages  and 

determine  tne  mass  of  the  ion  which  sucessfully  traverses  the  length 
of  the  instrument.  In  the  typical  operation  of  the  instrument  (1) 
and  are  allowed  to  vary  while  their  ratio  is  held  constant  and  tne 
number  of  ions  which  sucessfully  traverse  the  filter  are  recorded. 

The  mass  resolving  power  of  a  mass  filter  at  mass  m  is  defined  by 
m/am  where  Am  is  a  measure  of  the  smallest  difference  in  mass  at  mass  m 
that  can  be  detected  by  the  instrument  at  a  specified  valley  (e.g.  10% 
peak  height)  between  peaks  of  equal  heights.  The  transmission  of  a  mass 
filter  tuned  to  mass  m  is  defined  as  the  ratio  of  the  number  of  ions  with 
mass  m  injected  into  the  mass  filter  to  the  number  which  arrive  at  the 
detector.  It  would  be  desirable  for  a  QMF  to  have  high  resolving  power 
and  high  transmission.  Unfortunately,  these  are  conflicting  requirements 
and  the  improvement  in  one  of  these  parameters  results  in  a  degradation 
of  the  other  parameter. 
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Figure  1.  Electrode  structure  of  the  quad- 
rupole  and  its  relationship  with  the  cartesian 
coordinate  system.  The  coordinate  system  is 
oriented  as  shown  but  translated  so  that  the 
origin  is  centered  between  the  pole  pieces  as 
in  figure  3.  The  ions,  which  move  in  a  vacuum 
are  introduced  with  a  small  velocity  at  x=0, 
y=0  in  the  z  direction.  By  imposing  a  voltage 
difference  with  DC  and  AC  components  on  adjacent 
electrodes  the  device  is  tuned  so  that  only 
ions  with  the  proper  mass  to  charge  ratio 
traverse  the  length  of  the  filter  where  they 
are  detected.  Ions  with  a  different  mass  to 
charge  ratio  collide  with  the  electrode  walls 
and  are  lost. 

Figure  2.  Inverted  pendulum.  The  rigid  rod 
OP  of  length  Z  is  free  to  rotate  in  the  x-z 
plane.  The  z-axis  is  oriented  so  it  points 
away  from  the  center  of  the  earth  and  the 
supnort  point  0  has  an  acceleration  of 
-baj2  cos  tut.  For  a  high  enough  frequency  u) 
the  inverted  pendulum  has  the  remarkable 
property  of  being  stable  in  a  vertical 
position  with  its  center-of  mass  above  its 
support  point. 


The  research  reported  here  is  motivated  by  the  following  question: 

Can  the  shape  of  the  pole  pieces  or  the  time  dependence  of  the  applied 
potential  difference  be  changed  in  such  a  way  as  to  improve  the  resolution- 
transmission  relationship?  In  part  2  of  this  paper,  differential  equations 
of  ion  motion  are  derived  for  geometries  which  generalize  the  popular 
quadrupole  geometry.  By  integrating  these  equations  of  motion,  the 
suitability  of  these  geometries  for  mass  spectrometry  can  be  investigated. 
Should  calculation  indicate  a  particular  geometry  is  suitable,  the 
question  of  whether  or  not  it  yields  improved  performance  can  be  investi¬ 
gated  experimentally.  Part  3  of  this  paper  deals  with  the  relationship 
between  the  QMF  and  the  inverted  pendulum.  The  inverted  pendulum,  which 
consists  of  a  rigid  rod  free  to  rotate  in  a  vertical  plane  and  whose 
point-of-support  is  made  to  vibrate  vertically  has  the  remarkable  property 
of  being  stable  in  a  vertical  position  with  its  center  of  mass  above  its 
support  point,  providing  the  support  point  oscillates  above  a  critical 
frequency.  In  part  3  of  this  paper  it  is  shown  that  the  differential 
equations  for  an  ion  traversing  a  QMF  and  the  differential  equations 
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describing  the  motion  of  an  inverted  pendulum  can  be  cast  into  the  same 
form.  In  this  way  the  inverted  pendulum  is  shown  to  be  a  raecnanical  analog 
of  the  QMF.  This  aspect  of  the  work  is  significant  for  three  reasons: 

1)  The  motion  of  an  ion  as  it  traverses  a  QMF  cannot  be  physically  seen. 

By  observing  the  motion  of  an  inverted  pendulum  one  may  better  appreciate 

the  subtleties  involved  in  the  way  a  quadrupole  works  and  perhaps  obtain  a 
better  sense  of  what  is  happening  in  this  device.  This  may  be  of  value  to 
the  person  who  tries  to  design  a  better  mass  filter.  2)  The  inverted 
pendulum  has  pedagogical  value  for  explaining  the  operation  of  a  QMF  to 
students  and  has  already  been  used  for  this  purpose  by  one  of  the  authors 

(A.L.Y).  3)  By  showing  that  the  inverted  pendulum  and  QMF  are  related 

subjects,  it  is  possible  that  studies  on  one  of  them  may  shed  light  on 
the  other. 

The  inverted  pendulum  has  been  described  by  several  authors  (3-10)  . 
Corben  and  Stehle  (3)  have  modeled  the  inverted  pendulum  as  a  rigid  rod 
while  Landau  and  Lifshitz  (4)  have  modeled  it  as  a  point  mass  attached  to 
a  rigid  massless  rod.  Kapitza  (5)  has  discussed  the  case  of  the  pendulum 
when  the  amplitude  of  the  sinusoidal  oscillation  of  the  suspension  point 
is  small  in  comparison  to  the  length  of  the  pendulum.  He  has  also  reported 
what  appears  to  be  the  first  operative  pendulum  with  an  oscillating 
support  point  (normal  and  inverted).  Kalmus  (6)  has  discussed  the  theory 
of  the  inverted  pendulum  with  triangular  excitation  of  the  support  point, 
and  has  reported  an  operative  inverted  compound  pendulum.  Mitchell  (7) 
has  applied  a  method  of  averaging  to  find  the  staDility  of  an  inverted 
pendulum  for  small  amplitude  high  frequency  oscillation,  and  Howe  (8)  has 
described  a  theory  for  the  stability  of  an  inverted  pendulum  driven  by 
oscillations  of  a  random  nature.  A  square-wave  model  of  the  inverted 
pendulum  has  been  recently  described  by  Yorke  (9).  In  a  previous 
publication  CIO),  Friedman  et  al.  have  observed  that  the  inverted  pendulum 
is  a  mechanical  analog  of  the  QMF, 

Partial  derivations  of  the  differential  equations  for  ion  motion  in  a 
QMF  are  given  by  Campana  (1)  and  also  by  Dawson  (2).  A  derivation  (start¬ 
ing  from  the  Lorentz  force  law  and  the  most  general  form  of  the  Maxwell 
equations)  for  ion  motion  in  multipole  fields  with  a  discussion  of 
various  approximations  has  been  given  by  Friedman  et  al.  (11). 

2.  ION  MOTION  IN  MULTIPOLE  FIELDS 

a.  Fundamental  physics  of  ion  motion.  The  force  on  an  ion  with 
charge^e  moving  with  velocity^v  through  space  in  which  there  is  an  electric 
field  E  and  a  magnetic  field  B  is  given  by  the  Lorentz  force  law 

F  =  e(E  +  v  x  b)  .  (1) 

The  differential  equation  which  describes  the  motion  of  an  ion  through  a 
field  is  obtained  from  Newton's  law  of  motion 
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7  + 
d  r 
m  — 

dt 


(2) 


It  is  necessary  to  compute  the  E  and  B  fields  caused  by  the  time  varying 
potential  imposed  on  the  pole  pieces  to  find  the  motion  of  the  ions  from 
equations  (1)  and  (2).  These  are  computed  from  Maxwell's  equations.  If 
the  charges  on  the  ions  and  their  associated  current  densities  can  be 
neglected  for  the  purpose  of  computing  the  E  and  B  fields  acting  on  them, 
then  in  the  vacuum  region  where  the  ions  are  moving  Maxwell's  equations  are: 


-4-  -4> 

7*E  = 


0 


V-B  =  0 


V  x  E  =  - 


3B 

(EoE) 


(mj  ; 


If  the  E  and  B  fields  are  not  changing  with  time  then  Maxwell's  equations 
(Ml)  become: 

^•E  =0  ^  x  £  =  0 


^•B  =  0  V  x  B  =  0  (M2) 

The  applicability  of  the  two  forms  of  Maxwell's  equations  (Ml)  and 
(M2)  to  ion  motion  in  a  mass  filter  merits  discussion.  Maxwell's  equations 
(Ml)  are  a  good  approximation  in  the  vacuum  region  enclosed  by  the  pole 
pieces  of  a  mass  filter.  As  pointed  out  in  the  introduction,  the^voltages 
on  the  electrodes  change  at  radio-frequencies  resulting  in  B  and  E  fields 
which  also  change  at  the  same  frequencies.  Therefore,  it  is  not  obvious 
that  Maxwell's  equations  (M2)  can  be  used  for  the  purpose  of  computing 
ion  motion.  The  use  of  equations  (M2)  for  this  purpose  can  be  understood 
through  the  following  discussion.  Equations  (M2)  predict  (12)  that 
B0  =  E0/c  where  B0  and  E0  are  the  maximum  amplitudes  of  the  B  and  E  fields 
and  c  is  the  speed  of  light.  Thus  the  B  field  in  the  Lorentz  force  law 
(1)  exerts  a  negligible  force  on  the  ion  compared  to  the  E  field,  unless 
the  ion  is  moving  at  a  speed  comparable  to  the  speed  of  light.  Campana 
(1)  has  shown  that  the  fastest  ions  in  a  quadrupole  mass  filtar  are  moving 
slowly  (v/c  <  10  4)  compared  to  the  speed  of  light,  even  under  extreme^ 
operating  conditions.  Therefore,  it  is  not  necessary  to  consider  the  B 
fields  any  further  for  the  purpose  of  computing  the  force  on  an  ion.  The 
Lorentz  force  law  (1)  reduces  to  F=eE,  while  equations  (M2)  become: 


-4-  -4- 

V*E  =  0 


^  x  I  =  0 


(M3) 


The  second  of  these  equations  (M3)  is  valid  only  for  an  alternating 
potential  which  has  a  frequency  v  low  enough  such  that 
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X  =  -  »  l  (3) 

v 

That  is,  the  wavelenth  X  associated  with  the  electromagnetic  wave  must  be 
much  greater  than  the  length  t  of  the  electrode  structure.  The  necessity 
of  equation  (3)  can  be  understood  in  the  following  way.  If  the  wavelength 
does  not  satisfy  equation  (3)  there  is  a  possibility  of  standing  waves 
(just  as  in  the  case  of  sound  waves  in  an  organ  pipe)  and  these  are  not 
predicted  by  equations  (M3).  Typically,  the  length  of  the  mass  filter  is 
less  than  one  meter  and  the  frequencies  used  (1)  are  well  below  10  MHz. 
Therefore  equation  (3)  is  well  satisfied  and  equations  (M3)  can  be  used  to 
compute  the  ion  trajectory. 

An  E  field  which  satisfies  equation  (M3)  can  be  found  with  the  aid  of 
two  mathematical  theorems,  ^tokes  theorem,  valid  for  any  continuous, 
differentiable  vector  field  A,  asserts  that 

!  I c  (V  x  A)*dS  =  (j)  A*ar  (4) 

—h 

where  dr  is  a  vector  on  the  circumf err  nee  of  the  enclosed  area  S  and  whose 
direction  is  in  the  direction  of  integration.  The  vector  dS  has  a 
magnitude  equal  to  the  area  dS  and  its  direction  is  normal  to  dS.  The 
ambiguity  in  th<|  two  directions  normal  to  dS  is  resolved  by  the  right-hand 
rule,  that  is  dS  points  m  the  direction  of  the  thumb_^with  the  other 
fingers  of  the  right  hand  pointing  in  the  direction  dr.  In  equation  (4), 
the  integral  on  the  right  is  over  the  entire  circumference  and  the  integral 
on  the  left  is  over  the  area.S  (not  necessarily  plane)  enclosed  by  the 
circumference.  Stokes  theorem  enables  a  two-dimensional  integral  over  a 
surface  which  has  the  form  on  the  left  side  of  equation  (4)  to  be 
evaluated  by  the  simpler  one-dimensional  integral  on  the  right-hand  side  of 
equation  (4).  Application  of  Stokes  theorem  to  the  second  of  equations 
(M3)  implies 

&  E*dr  =  0  (5) 

over  any  closed  path  which  defines  a  surface  over  which  ^  *  E  =  0.  The 
second  mathematical  theorem  asserts  that  any  vector  field  E  which 
satisfies  equation  (5)  around  every  closed  path  can  be  represented  as  the 
divergence  of  a  scalar  field  4>(x,y,z),  that  is 

E  =  -tf<J)(x,y,z)  •  (6) 

Substituting  E  from  equation  (6)  into  the  first  of  equations  (M3)  gives 
V*V<Kx,y,z)  =  V2t|)(x,y,z)  =  0  .  (7) 

The  prescription  for  finding  the  differential  equations  of  ion  motion 
for  any  configuration  of  pole  pieces  can  now  be  given.  Let  f  (x,y,z) 
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denote  the  equation  for  the  surface  S.  of  conductor  i.  The  differential 
equation  of  motion  can  be  found  by  first  obtaining  a  potential  function 
<j)(x,y,z)  which  satisfies  Laplace's  equation  (7)  and  which  also  satisfies 
the  boundary  conditions 


<f>(x,y,z) 


(8) 


•  1 

Generally,  it  is  difficult  to  find  a  function  which  satisfies 
equations  (7)  and  (8).  The  function  <J>(x,y,z)  can  sometimes  be  found 
analytically  using  the  method  of  images  (12-14),  or  by  complex  variable 
theory  (12-15),  or  by  infinite  series  expansions  (12-14,16)  or  it  can  be 
evaluated  numerically  using  digital  computers  (17-19)  .  Once  the  potential 
<f>(x,y,z)  is  known,  the  differential  equations  of  motion  can  be  found 
using  Newton's  law 


F  =  m  — -  =  -etfip  (x ,  y ,  z) 
dt“ 


(9) 


As  discussed  earlier,  this  method  is  approximately  valid  for  time-varying 
potentials 


h  ‘  V  ‘  *AC  C°S  Ut  (10) 

providing  the  potential  does  not  change  too  rapidly  with  time,  i.e.  subject 
to  the  limitations  of  equation  (3).  Physically,  equation  (10) 
corresponds  to  subjecting  the  electrodes  to  a  voltage  which  has  DC  and  AC 
components. 

The  remainder  of  this  paper  will  consider  only  infinitely  long 
conductors  whose  shape  does  rot  vary  with  z.  Although  real  electrode 
structures  must  necessarily  be  of  finite  length,  this  approximation  is 
made  because  it  captures  the  essential  features  for  the  operation  of 
electric  RF  devices  while  avoiding  some  of  the  mathematical  complexities. 
With  this  approximation  the  potential  for  such  a  conductor  configuration 
does  not  depend  on  z  i.e,  ^(x.y.z)  =  <J)(x,y)  and  hence  the  force  in  the 
z=direction  is  zero  i.e.  m  z  =  0.  The  equation  of  motion  in  the 
z-airection  can  be  integrated  directly 

z  =  v  t  +  z0  .  (11) 

z0 

Here  v  and  zD  are  the  z-component  of  velocity  and  the  z-coordinate 
respectively  at  t=0. 

b.  Method  of  obtaining  potential  function  and  equations  of  ion 
motion  in  multipole  fields.  The  approach  of  the  previous  section  was  to 
find  the  potential  function  <})(x,y,z)  for  a  specified  set  of  conductors  at 
potentials  <j>  .  An  alternative  approach  is  to  find  a  function  Un(x,y) 
which  satisfies  Laplace's  equation  (7)  and  from  this  function  determine  the 
equipotential  surfaces  of  the  pole  pieces.  If  conductors  with  given 
applied  potentials  are  fashioned  which  coincide  with  the  equipotential 
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surfaces,  a  potential  function 

4>(x,y)  =  K  U  (x, y)  (12) 

n  n 

can  be  found  which  satisfies  Laplace's  equation  and  the  boundary 
conditions  (8)  by  simply  choosing  the  constant  K  properly.  Physically,  K 
is  related  to  the  dimension  of  the  electrode  structure  and  the  applied 
electrode  potentials;  it  is  chosen  such  that  c})(x,y)  matches  the  known 
potentials  at  the  equipotential  surfaces. 

Functions  which  satisfy  Laplace's  equation  may  be  found  by 
application  of  the  theory  of  complex  variables.  If  the  complex  variable 
z  =  x  +  iy  is  raised  to  an  integral  power  of  n,  then  the  result  can  be 
expressed  as  the  sum  of  two  functions,  a  real  U  (x,y)  and  and  imaginary 
vn(x,y)  part: 

(x  +  iy)n  -  U  (x,y)  +  iV  (x,y)  .  (13) 

n  n 

The  integer  n  defines  the  order  of  the  multipole  field.  From  complex 
variable  theory  z  is  known  to  be  analytic  and  so  the  function  U  (x,y) 
satisfies  (15)  Laplace's  equation  (7). 

The  method  described  here  for  obtaining  the  potential  function  and  the 
differential  equation  of  motion  can  be  summarised  as  follows:  raise  the 
complex  variable  z  to  an  integer  power  n  to  find  a  potential  function  which 
satisfies  Laplace's  equation  (7)  and  the  boundary  conditions  (8);  then  ion 
motion  is  found  from  Newton's  Law  of  motion  (9). 

c.  Application  to  some  multipole  geometries.  The  method  described  in 
section  2b  will  now  be  applied  to  the  hexapole  geometry.  The  function  z 
is  evaluated  for  n=3  to  obtain 

U3(x,y)  =  x3  -  3xy2  . 

2 

Since  7“U^(x,y)  =  0,  the  function  is  the  basis  for  a  possible 

potential  function  in  charge  free  space.  Its  equipotential  surfaces  are 
illustrated  in  figure  4.  Each  of  the  pole  pieces  have  the  same  shape  and 
so  the  structure  is  unchanged  under  a  60  degree  rotation.  The  pole  piece 
lying  on  the  positive  x-axis  is  chosen  to  have  positive  applied  potential 
-KJ> o / 2  to  comply  with  convention  (2)  and  the  adjacent  electrodes  are 
chosen  to  have  potential  with  the  same  magnitude  but  opposite  sign.  The 
potential  function  d)(x,y)  which  satisfies  Laplace's  equation  and  the 
boundary  conditions  is 

<f>(x,y)  =  (x3  -  3xy2)  (16) 

2r, 

as  is  readily  verified.  With  <j>0  given  by  equation  (10),  Newton's  law 
becomes : 


7 


i.  iXi-'Jl.n-uit'i  ,  v—  aU  1  i-  «  X  i— iXWsJI  i-i  X 


FIGURE  3.  QUADRUPOLE  GEOMtTRY. 
EQUATION  OF  THE  POLE  PIECES; 

<p(x,  y)  =  ((fso/lrltx  -y1). 


FIGURE  5,  OCTAPOLE  GEOMETRY. 
EQUATIONS  OF  THE  POLE  PIECES: 

<p(x,  y)  =  (<po/2roXx*  -  6  x2y2  +  /). 


FIGURE  4.  HEXAPOLE  GEOMETRY. 
EQUATIONS  OF  THE  POLE  PIECES: 

tp(x,  y)  —  (©o/2ro)U3  -  3xy2). 


FIGURE  6.  DECAPOLE  GEOMETRY. 
EQUATIONS  OF  THE  POLE  PIECES: 

<p(x,y)  =  ((p o/2roXxs  —  lOxV  +  5-xy*). 
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The  method  just  illustrated  for  n=3  determined  the  shape  of  the  pole  pieces 
(equation  (16))  and  the  equations  of  ion  motion  (equation  (17))  for  the 
hexapole  geometry.  A  similar  procedure  can  be  used  for  n  any  integer. 

The  geometry  of  the  pole  pieces  when  n=2  is  the  well-known  quadrupole 
geometry  illustrated  in  figures  1  and  3.  The  differential  equation  for  ion 
motion  in  this  case  is: 
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cos  CJt)  x  =  0 


cos  tut)  y  =  0 
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The  geometry  of  the  pole  pieces  when  n=4  is 
The  differential  equation  for  ion  motion  in  this 

d2x  ,  2e  .  .  ,  3  2. 

m  — 2  +  ~b  ^DC  "  ^AC  C°S  (jQC^  ^X  “  3xy  ^  = 

dt  rQ 

2 

m  _  2e_  ^  cog  ujt)  (3x2y  _  y3^  = 

dt  r0 


illustrated  in  figure  5. 
case  is: 


0 
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The  geometry  of  the  pole  pieces  when  n=5  is  illustrated  in  figure  6. 
The  differential  equation  for  ion  ion  motion  in  this  case  is: 


d2x  5e 

"  7J  +  ~5  WDC 
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<f)  cos  out)  (x  -6xy  +  y)=0 


3  3 

p  cos  OJt)  (x  y  -  xy  )  =  0 

A.  Li 


(21) 


The  method  illustrated  in  these  examples  can  easily  be  used  for  finding  the 
differential  equations  of  ion  motion  for  n>5. 


3.  THE  INVERTED  PENDULUM  AND  THE  QUADRUPOLE  MASS  FILTER 

a.  Theory  of  the  inverted  pendulum.  In  this  section  the  differential 
equation  of  a  rigid  pendulum  with  vibrating  support  point  is  derived. 

Figure  2  shows  the  pendulum  which  is  a  thin  solid  rod  of  length  t  and  mass 
m  tilted  at  angle  0  with  respect  to  the  positive  z-ax.is  and  whose  point  of 
support  is  moving  with  speed  Z  in  the  positive  z  direction.  The  equation 
of  motion  for  the  pendulum  is  found  by  use  of  the  Lagrangian  function  L 
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where  L(q^_,q^)  =  ^(qk,qk)  ~  V(q  )  and  q^,  q  ,  T  and  V  are  the  generalized 
coordinaCes,  velocities,  non-relativistic  kinetic  and  potential  energies 
respectively.  The  differential  equation  of  motion  is  given  by  Lagrange's 
equation 


3L(V\> 


3L(Vqk> 


Consider  the  motion  of  the  k  mass  point 


at  distance  /  from  the 
m,  is  the  vector  sum  of 


support  point.  The  velocity  v^  of  the  mass  point  is  the  vector  sum  of 
the  velocity  due  to  rotation  and  the  velocity  due  to  motion  of  the  support 


point 


-  _  ^k  de_  ?  dzk  de  •  r 

Vk  d0  dt  1  (d0  dt  z)  k 


The  geometrical  relationships 


xk  =  ek  sin 


z,  =  L  cos  9 
k  k 


allow  v^  to  be  expressed  in  terms  of  the  generalized  coordinate  9  and 
generalized  velocity  0 


vk  =  £k  0  cos  9  i  +  (Z  -  £k  6  sin  0)  k  .  (23) 

Here  i  and  k  are  unit  vectors  along  the  x  and  z  axis  respectively. 
Equation  (23)  and  the  definition  of  kinetic  energy, 

T  '  2  (  Vk  ’  1  l  \  Vic  • 

allow  the  kinetic  energy  to  be  expressed  in  terms  of  the  generalized 
coordinates  and  velocities 

T  =  \  (l  \-ek)  02  -  (l  m k£k)  0  Z  sin  9  +  \  m  Z2  .  (24) 

k  k 

As  seen  from  figure  2,  the  potential  energy  V  of  the  pendulum  in  terras  of 
each  mass  point  is  a  decreasing  function  of  9  given  by 

v  =  (.1  8  cos  9  (25) 

k 

where  g  is  the  acceleration  of  gravity.  Equations  (24)  and  (25)  may  be 
simplified  by  using  the  definitions  for  the  moment  of  inertia  I  of  a  body 
rotating  about  the  support  point  and  the  coordinate  of  the  center  of 
mass  measured  from  the  support  point  along  the  axis  of  the  rod: 


=  1  l 


—  I  m  l 
ra  ,  k  k 
k 
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If  the  mass  is  distributed  symmetrically  about  its  mi^-point,  then  t  -ill 
and  Lagrange's  equation  can  be  written 

1  •?  1  .?  9  .  .  } 

L  =  -  I  f  +  y  m  Z"  -  m  y  0  Z  sin  0  -  mg  —  cos  0  .  (26) 

With  q=0 ,  Lagrange's  equation  (22)  applied  to  equation  (26)  yields  the 

differential  equation 

l  "  f 

10  -  m  —  Z  sin  0  =  mg  —  sin  0  .  (27) 

Two  simplifications  of  equation  (27)  can  be  made.  First,  the  moment 
of  inertia  for  a  uniform  thin  solid  rod  is  easily  computed  (10) 

I  £2 

Second,  the  approximation  sin  0  =  0  is  made  since  it  is  anticipated  that 
the  oscillations  of  the  pendulum  will  be  stable  near  0=0  and  so  sin  0 
will  be  small.  Making  these  two  substitutions,  equation  (27)  becomes 


0  -  fr  (g  +  2)  0  =  0  . 


Equation  (28),  the  principal  result  of  section  (3c),  is  seen  to  be 
reasonable  if  the  cases  Z=0,  Z=-g  and  Z=+g  are  considered.  For  the  case 
Z=0,  equation  (28)  asserts  that  0  will  grow  exponentially  with  time 
constant  t  =  ( (2£) / (3g) ) *.  This  is  in  agreement  with  the  intuitive 
expectation  that  an  inverted  pendulum  without  oscillations  is  unstable. 

For  the  case  Z=-g,  the  support  point  is  moving  down  with  a  constant 
acceleration  g.  Equation  (28)  predicts  that  0=0  and  this  is  in  agreement 
with  intuitive  expectations  since  for  a  freely  falling  pendulum  the 
effective  forces  acting  to  change  the  variable  0  are  expected  to  vanish. 

For  the  case  Z=+g  the  pendulum  support  point  is  accelerating  upward  with 
constant  acceleration  g.  Equation  (28)  predicts  that  the  motion  for  the 
variable  0  is  the  same  as  that  of  a  stationary  pendulum  sitting  in  a 
gravitational  field  of  strength  2g  and  this  too  is  in  agreement  with 
intuitive  expectations. 

b.  The  Mathieu  equation.  Subsequently  it  will  be  seen  that  the 
equation  describing  ion  motion  in  a  quadrupole  and  the  equation  describing 
the  inverted  pendulum  can  both  be  put  in  the  cannonical  form  of  the 
Mathieu  equation.  The  purpose  of  this  section  is  to  describe  this  equation 
and  its  stability  characteristics. 

The  canonical  form  of  the  Mathieu  equation  is 

,2 

—  +  (a  -  2q  cos  2£)u  =  0  .  (29) 

d  C 

Here  a  and  q  are  real  constants,  positive  or  negative,  which  are  not 
functions  of  the  dependent  variable  u  or  the  independent  variable  £j. 


11 


w 


■  Has 
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i  /  Figure  7.  The  a-q  stability  diagram. 

Si  ./  Equation  (29)  defines  the  parameters  a  and 

q  which  appear  in  the  Mathieu  equation, 
a  For  a  and  9  values  which  fall  in  the 

k  n.  shaded  area  the  value  of  u  in  equation  29 

! *->,  j  V  jq  ^  *  remains  bounded  as  z,  increases  independ- 

\  ent  of  the  initial  conditions.  Other 

\  \  values  of  a  and  q  result  in  unbounded 

.5!  solutions. 

Although  the  factor  -2  multiplying  q  and  the  2  in  the  argument  of  the 
cosine  function  might  look  strange,  this  Is  the  standard  for  of  this^^  ^3 
differential  equation  whose  properties  have  been  extensively  studied.  ’ 

The  solution  of  equation  (29)  will  be  a  function  of  the  form 

u  =  u(£;  a,q,u0,u0)  •  (3°) 

Equation  (30)  expresses  mathematically  the  observation  that  the  function 
u  depends  on  the  independent  variable  £,  on  the  parameters  a  and  q,  on  the 
initial  value  of  u  denoted  by  u0.  and  on  the  initial  value  of  du/d£ 
denoted  by  u0.  The  function  (30) # evaluated  for  a  particular  value  of  a 
and  q  is  a  function  of  q,  u0  and  uQ.  The  result  of  detailed  mathematical 
analysis  (22)  is  that,  depending  on  the  values  of  the  parameters  a  and  q, 
the  function  u  becomes  unbounded  or  remains  bounded  as  c,  increases  and 
that  these  values  of  a  and  q  are  independent  of  the  initial  conditions 
u0  and  uQ.  Figure  7  shows  a  portion  of  a  stability  diagram  in  a-q  space. 
The  shaded  area  results  in  "stable"  trajectories  by  which  it  is  meant  that 
u  remains  within  finite  bounds  as  E,  increases  independent  of  initial 
conditions . 

c.  Analogy  between  inverted  pendulum  and  quadrupole  mass  f ilter .  In 
this  section  it  is  shown~that  both  the  inverted  pendulum  and  the  QMF  can 
be  expressed  in  the  cannonical  i  jrm  of  the  Mathieu  equation  cud  the 
analogs  will  be  identified. 

In  equation  (28)  Z  is  the  acceleration  of  the  support  point.  If  the 
support  point  is  driven  sinusoidally  with  frequency  w,  then  the  origin  of 
time  can  be  chosen  so  that  Z  =  b  cos  uJt,  In  that  case,  equation  (28) 
becomes 

2 

(g  -  b  cos  ujt)  8  =  0  .  (31) 

dt“  2 

Define  four  dimensionless  parameters 


u  =  9 

(b) 

a  =  -  ^ 

£u) 

11 

IS 

(d) 

3b 

q  =  -  - 
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and  use  these  to  replace  the  variables  6  and  t  in  equation  (31)  to  get  the 
cannonical  form  of  the  Mathieu  equation  (29)  . 

Note  that  equations  (19a)  and  (19b)  are  essentially  the  same  aside 
from  a  change  in  sign  and  so  without  loss  of  generality  only  equation 
(19b)  is  considered  in  drawing  the  analogy  with  the  inverted  pendulum. 
Define  four  dimensionless  parameters 


u  =  y 


2  2 
mu)  r0 


(c)  £  = 


q  =  - 


2  2 
mu)  rD 


and  use  these  to  replace  y  and  t  in  equation  (19b)  to  get  the  cannonical 
form  of  the  Mathieu  equation  (29). 

From  equations  (32)  and  (33)  the  analogous  quantities  can  be  ident¬ 
ified.  Multiplying  numerator  and  denominator  of  equation  (32b)  by  m,  the 
mass  of  the  pendulum,  and  comparing  equations  (32 d)  and  (33b),  aside  from 
some  unimportant  numerical  factors,  the  quantity  mg  of  the  pendulum  is 
analogous  to  (“'('qq) /r°  of  Che  QMF ;  ml  of  the  penaolum  Lj  analogous  to  mrQ 
of  the  QMF  and  w  of  the  pendulum  is  analogous  to  oj  ui  the  QMF. 

Multiplying  numerator  and  denominator  of  equation  (32d)  by  mur  and 
comparing  equations  (32d)  and  (33d),  aside  from  some  unimportant 
numerical  factors,  the  quantity  bmui2  of  the  pendulum  is  analogous  to 
(e^Ac) /c j  of  the  QMF;  ml  of  the  pendulum  is  analogous  to  mr0  of  the  QMF 
and  to  of  the  pendulum  is  analogous  to  u)  of  the  QMF.  In  these  analogies, 
all  the  quantities  associated  with  the  pendulum  have  the  same  physical 
dimensions  as  the  corresponding  quantities  associated  with  the  QMF.  Also, 
the  RF  source  of  the  QMF  and  the  oscillating  source  for  the  support  point 
in  the  inverted  pendulum  perform  analogous  functions  of  producing  a 
stability  mechanisim  for  each  system.  The  destabilizing  terms  of  each 
system  are  g,  the  acceleration  of  gravity  for  the  inverted  pendulum, 
and  ,  the  applied  DC  voltage  of  the  QMF,  The  quantities  which  counter 
these  destabilizing  terms  are  bio2  ,  the  acceleration  amplitude  for  the 
support  point  of  the  inverted  pendulum  and  4>  ,  the  RF  amplitude  of 

the  QMF.  L 

4.  CONCLUSIONS. 

Explicit  differential  equations  in  closed  form  have  been  found  for 
pole  piece  geometries  which  are  generalizations  of  those  used  in  the 
QMF.  For  arbitrarily  shaped  pole  pieces,  the  solution  of  Laplace's 
equation  cannot  be  written  in  closed  form  and,  thus  in  general  the  ion 
differential  equations  of  motion  cannot  b  written  in  closed  form.  The 
geometries  shown  in  figures  4,  5  and  6  are  analytically  useful  in 
determining  if  other  geometries  might  yield  mass  filters  with  improved 
resolution/ transmission  characteristics.  Computer  simulation  of  ion 
motion  in  these  geometries  is  a  subject  for  future  research. 


FRIEDMAN,  CAMP  ANA  &  YERGEY 


Although  not  reported  in  this  paper,  an  inverted  pendulum  has  been 
built  (10)  and  it  has  been  useful  in  emphasizing  to  students  the  non- 
intuitive  nature  of  the  mass  filtering  action  that  takes  place  in  a  QMF. 

In  reference  (10) ,  it  is  shown  conceptually  that  the  inverted  pendulum 
can  be  used  to  experimentally  measure  the  boundaries  between  the  stable 
and  unstable  regions  in  the  a-q  stability  diagram.  As  a  result,  the 
inverted  pendulum  aids  in  understanding  these  diagrams  which  are  crucial 
to  understanding  the  operation  of  a  QMF. 

The  principal  new  results  of  our  research  are:  1)  new  derivations  of 
the  equations  which  describe  ion  motion  in  a  QMF,  2)  the  equations  of  pole 
pieces  for  hexapole,  octupole  and  decapole  geometries  illustrated  in 
figures  4,  5  and  6,  3)  closed  form  differential  equations  of  ion  motion  in 
these  geometries  (equations  (17),  (20)  and  (21))  and  4)  the  realization 
that  the  inverted  pendulum  is  a  mechanical  analog  of  the  QMF.  _ 

Today  the  quadrupole  mass  analyser  is  the  most  widely  used  mass 
analyser  for  low  resolution  applications.  Typically,  mass  filters  are 
used  in  laboratory  environments  where  simplicity,  compactness,  economy, 
absence  of  magnetic  fields,  and/or  the  capability  of  fast  scan  rates 
(especially  for  chromatographic  combinations)  are  important  (1).  Thus, 
an  improvement  in  resolution/ transmission  qualities  of  electric  RF  mass 
filters  would  enable  this  widely  used  mass  spectrometer  to  take  a  larger 
role  in  analytical  laboratories  and  in  basic  research.  Other  design 
considerations  besides  resolution/ transmission  characteristics  are 
important.  If  mass  filters  could  be  made  which  are  smaller,  more  rugged, 
less  susceptible  to  vibration  with  reduced  power  requirements,  and  more 
easily  maintained,  there  might  be  greater  application  outside  the 
laboratory.  Possible  applications  include:  surveying  for  natural 
resources,  pin-pointing  sources  of  industrial  pollutants,  and  the 
detection  of  hazardous  materials  such  as  chemical  agents  and  explosives. 
Eventually  our  work  may  be  applied  to  the  solution  of  these  problems. 
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